Ultrafast electro-nuclear dynamics of H2 double ionization 
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The ultrafast electronic and nuclear dynamics of H 2 laser-induced double ionization is studied us- 
ing a time-dependent wave packet approach that goes beyond the fixed nuclei approximation. The 
double ionization pathways are analyzed by following the evolution of the total wave function during 
and after the pulse. The rescattering of the first ionized electron produces a coherent superposition 
of excited molecular states which presents a pronounced transient H + H~ character. This attosecond 
excitation is followed by field-induced double ionization and by the formation of short-lived autoion- 
izing states which decay via double ionization. These two double ionization mechanisms may be 
identified by their signature imprinted in the kinetic-energy distribution of the ejected protons. 

PACS numbers: 33.80.Rv, 33. 80. Eh, 42.50.Hz 
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Laser-induced non-sequential double ionization of 
atoms has been extensively studied since its discovery in 
the early 1980s [1]. This strong field effect as well as high- 
order harmonic generation (HHG) in rare gases Q gave 
birth to the so-called rescattering model which gives a 
unified picture of the atomic response 0]. Moreover, this 
model constitutes the natural framework for the produc- 
tion of attosecond pulses and for the control of ioniza- 
tion in atoms and molecules On the molecular side, 
new effects in strong laser fields involve the coupling of 
the electronic and nuclear motions such as bond soften- 
ing [5] or charge resonance enhanced ionization [6] . Re- 
cently, this coupling was shown to play a determinant 
role in HHG [7, 8, 9]. In particular, Lein predicted that 
more intense harmonics are generated in heavier molec- 
ular isotopes [7] . This effect has been demonstrated ex- 
perimentally using H2, D2, CH4 and CD4 [8]. In addi- 
tion, molecular vibrations were recently detected with a 
mA sensitivity using HHG in SFq [9|. Following a recent 
theoretical proposal , ground state vibrational wave 
packets of D2 where also mapped with 20 mA sensitivity 
using selective depletion via tunnel ionization . 

These recent advances show that nuclear motions are 
fast enough to exercise a decisive influence on the over- 
all electronic response of the molecule in strong fields 
although the vibrational time scale only belongs to the 
femtosecond domain. Moreover, the hydrogen molecule 
exhibits the shortest vibrational period: 7.5 fs. From 
the theoretical point of view it is therefore highly desir- 
able to develop a unified and realistic theoretical frame- 
work where the nuclear dynamics is treated at the same 
quantum level as the electronic excitation. By treating 
the internuclear distance R as a full quantum variable, 
we show that impulsive vibrational excitation also takes 
place in H2. In addition, the rescattering dynamics is 
basically a two-electron problem which demands a two- 
active-electron approach. This approach allows for the 



description of the interplay between electronic rescatter- 
ing and nuclear vibrations. Using a single-optical-cycle 
pulse, we show that the electron rescattering leads dur- 
ing the pulse to field-induced double ionization, and after 
the end of the pulse to auto-double ionization. 

Using the splitting technique of the short-time prop- 
agator, we solve the time-dependent Schrodinger equa- 
tion for H2 in the linearly polarized classical electric field 
E(t) = Eq fit) cos(cot + if) e, where uo denotes the angu- 
lar frequency, Eq the field amplitude, and ip the carrier- 
envelope offset phase. The pulse shape, of total dura- 
tion r, is f(t) = sin 2 (7rt/r). The frequency uo corre- 
sponds to a central wavelength of 800 nm. The inter- 
nuclear coordinate R is constrained along the field po- 
larization vector e and the two electrons, of coordinates 
7*1 — — z>i 6, are assumed to oscillate along the same axis. 
This one-dimensional approach was shown recently to 
represent faithfully most of the electronic dynamics [12| . 
The wave function ^(R, z±, Z2, t) of H2 is taken initially 
as the ground state of the soft- Coulomb potential 
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which takes into account the electrostatic repulsion and 
attraction between all particles [13[ . In order to mimic 
the dynamics of the real H2 molecule, we have intro- 
duced two softening parameters, a(R) and (3(R), which 
vary slowly with the internuclear distance. These param- 
eters have been adjusted at each internuclear distance to 
reproduce accurately the potential curves of H2 and H^ 
ground electronic states. All spectroscopic constants and 
ionization potentials of our model H2 and Hj molecules 
therefore match very closely the tabulated ones. 

During the pulse, the Hamiltonian includes the radia- 
tive coupling [z\ + Z2)E{t) in the dipole length represen- 
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FIG. 1: (Color online) (a) Theoretical and experimental pro- 
ton kinetic energy spectra at 4.5 x 10 14 W/cm 2 . The mea- 
sured pulse duration (FWHM) is 10±2fs and the calcula- 
tion is for 11.7 fs. (b) Theoretical spectrum at 5xl0 14 W/cm 2 
and for r — 2.67 fs. The dotted and dashed lines correspond 
to recollision-induced field- assisted double ionization (RFDI) 
and to recollision-induced auto-double ionization (RADI). 



tation. The outgoing ionization flux is summed to ex- 
tract the single and double ionization probabilities. For 
this purpose, we partition the plane (zi,z 2 ) in three re- 
gions To, Ti and T 2 corresponding respectively to H2, 
Hj, and H^. Double ionization occurs in the asymp- 
totic region T2 = {|^i|, \z 2 \ > za] with za = 20 ao ^> R. 
The single ionization region Ti is defined as {\zi\ < za, 
\zj\ > ^a}, and the neutral H 2 molecule is found in the 
region Tq = {|^i|,|^2| ^ ^a}- The total wave function 
zi, z 2 , t) is propagated on a spatial grid of (512) 3 
points with z max = 100 ao, R max = 10 ao and 5t = 1 as. A 
simple mapping relates the proton kinetic energy distri- 
bution S(E,t) to P 2 (R,t) = JJ T2 \^(R,z 1 ,z 2 ,t)\ 2 dz 1 dz 2 
using E = 1/ (2R) and the requirement of particle con- 
servation P 2 (R,t) dR = S(E,t) dE. This energy distri- 
bution is accumulated over the entire time propagation 
to obtain the kinetic energy release spectrum S(E) which 
is measured in experiments. Fig.[TJa) shows that a good 
agreement is obtained for the Coulomb peak in the en- 
ergy range 3 eV ^ E ^ 8 eV between this model and the 
experiment The low energy peak (E ^ 3eV), which 
results from the photodissociation of in H + H + , is 
not reproduced by this Coulomb mapping. 

A simplified analysis of the various double ioniza- 
tion mechanisms can be performed using the shortest 
pulse containing a complete oscillation of the field, with 
r = 2tt/uj = 2.67 fs and tp = —tt/2. The associated elec- 
tric field E(t), which is shown Fig.[2fa) as a thin dotted 
curve, presents two symmetric maxima pointing in op- 
posite directions. The time variations of the single (SI) 
and double (DI) ionization probabilities are also shown 
Fig.^a) for a peak intensity of 5xl0 14 W/cm 2 . Single 
and double ionization take place by successive bursts. 
The probability of single ionization is simply rising just 
after each maximum of the field, and the ionized electron 
reaches the asymptotic region in At ~ 700 as. The evolu- 
tion of the double ionization probability is more complex. 
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FIG. 2: (Color online) (a): Single (SI) and double (DI) ion- 
ization probabilities as a function of time at 5xl0 14 W/cm 2 
and for r — 2.67 fs. The electric field E(t) is superim- 
posed on this graph. The vertical line marks the end of 
the pulse, (b), (c) and (d): Nuclear probability distribu- 
tions Pk(Rjt) — ff rk \ 1 $/(R 1 zi 1 Z2,t)\ 2 dzidz2 as a function 
of internuclear distance R and time t in the three regions Fk 
corresponding to H2 (b), (c) and (d). 



Indeed, double ionization takes place in four successive 
bursts which differ considerably by their magnitudes and 
time scales. Since the first appearance of double ioniza- 
tion (label (D, Fig. 0(a)) follows the first maximum of the 
field by about 700 as, it is clearly due to laser-induced di- 
rect double ionization. The second burst (©) just follows 
the second maximum of the field and thus probably arises 
from single ionization of the molecular ions already 
produced. Finally, we will show in the following that 
the two last and dominant pathways to double ionization 
(®,®) are induced by recollision and can be attributed 
to laser-assisted and to autoionization processes. 

Fig.[2][b), which represents the evolution of the wave 
packet in the To region, shows that, even for a pulse 
of total duration r = 2.67 fs, a significant nuclear mo- 
tion takes place in H2. At about 1.2 fs, the neutral 
molecule stretches from the equilibrium internuclear dis- 
tance 1.4 ao to 1.5 ao- This vibrational wave packet, co- 
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herent superposition of a few vibrational levels, arises 
from a small and early lochfrafi effect [10[ which is ampli- 
fied by a subsequent bond-softening like dynamics. The 
competition between these two mechanisms was moni- 
tored very recently in D2 [Tlj . 

The nuclear dynamics taking place during single and 
double ionization is even more significant. Figures[2jc) 
and [2](d) show that the chemical bond stretches up to 
1.7 do when the molecule is singly ionized and up to 3ao 
during double ionization. The first three bursts of double 
ionization are identified in Fig.[2](d) by the labels CD, ® 
and ®. These three processes take place at similar in- 
ter nuclear distances, 2 — 2.2 ao, close to the minimum of 
the potential associated to the ground electronic state of 
(2ao). The last mechanism (label ®) is taking place 
much later, t ~ 3 — 4.5 fs, after the end of the pulse, and 
around 3ao This distance is far beyond the equilibrium 
distances of H2 and Hj. In this reaction path, a dissoci- 
ation process is therefore already taking place during the 
second stage of ionization, before Coulomb explosion. 

A clear identification of the four mechanisms giving 
rise to Coulomb explosion of H2 is obtained by analyzing 
the electronic dynamics shown in Fig. [3] at specific times 
and internuclear distances. The times t and distances 
R we have chosen correspond to the four labels given in 
Fig.[2|(d). In agreement with the previous discussion, we 
observe that laser-induced direct double ionization takes 
place around 1.2 fs (white arrows in Fig.[3^a)). Indeed, 
the two electrons escape simultaneously from the nuclear 
attraction in the same direction (z\ and Z2 < 0), op- 
posite to the direction of the electric field. In addition, 
the electronic repulsion quickly washes out the electronic 
density around z\ — z^. The curvature of the two path- 
ways seen in this figure indicates that an energy exchange 
is taking place, inducing an acceleration of one of the 
electrons to the detriment of the other. Similar effects 
were seen earlier in the context of multi-electron ioniza- 
tion of atoms [15]. The electronic probability distribution 
of Fig. 02(b) also confirms that sequential double ioniza- 
tion occurs around 1.8 fs. In this case, the first electron 
has already left the molecular ion in the direction z < 
when the second electron is ejected. The driving electric 
field has then reversed sign, and this second electron is 
leaving the nuclei in the opposite direction. 

The nature of the two dominant pathways to double 
ionization is finally enlightened by the snapshots of the 
electronic density given in the last two panels of Fig. 03 
When the electric field changes sign (t > 1.33 fs), the 
part of the electronic wave function which is located in 
the single ionization region Ti, and which therefore rep- 
resents an ionized electron far from the bound electron 
of Hj, starts to be driven back towards the nuclei. This 
recollision process induces an interference between the 
incoming (recollision) and outgoing (ionization) electron 
waves. This interference is easily seen in the oscillatory 
behavior of the electronic density in Fig. 03(c). The two 
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FIG. 3: (Color online) Electronic probability distribution 
\^(R, zi, Z2,t)\ 2 (log scale) at different times and internuclear 
distances. The white arrows represent the photoelectron flux. 
The panels a-d correspond to the points 1-4 in Fig.[2jd). 



peaks identified in this figure, which are separated by 
Az ~ 5ao, will recollide successively with the electron 
that remained near the nuclei in a time interval of about 
115 as. Because of this interference pattern, the recol- 
lision takes place as a very fast succession of periodic 
electron impacts. Each impact induces a very specific 
electronic excitation in the central region where \z±\ and 
\z2 1 ^ R. Indeed, in the vicinity of the two nuclei the elec- 
tronic wave function does not show anymore the usual 
double peak structure of Eb's ground electronic state. 
This structure is now replaced by a dynamical corre- 
lated motion of the two electrons around the four crosses 
shown in this figure. The two signs correspond to 
the two attractive combinations where each electron is 
located on a different nucleus. The two 'x' signs corre- 
spond to the two ionic combinations H + H _ and H _ H + . 
Around the internuclear distance 2.2 ao, the energy gap 
between the ionic state and the ground electronic state 
of H 2 can be estimated from the energy difference be- 
tween the B 1 S+ and the X 1 S+ potential curves of H 2 
as A ~ 0.37 Eh. Considering the classical ponderomotive 
energy U p = e 2 E^^/Amuo 2 ~ 0.46 E^ for the present 
maximum value of the field E max ~ 0.08 E^/edo, the 
H + H _ spectral range is clearly accessible following the 
recollision process. In addition, recent studies [l6[ have 
shown that one of these ionic states is stabilized by strong 
fields, such that the ionic and the ground electronic states 
of H2 may experience an avoided crossing in an adiabatic 
picture. As a consequence, the ionic configuration H + H~ 
becomes the instantaneous ground state of the system 
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for large field amplitudes. We find here that this effect, 
which was first demonstrated in a two-electron model of 
H2 with frozen nuclei [l6[ , survives in a more realistic ap- 
proach which takes into account the dissociative dynam- 
ics of the molecule. For large field values, and because 
of the recollision, the two electrons are thus forced to lo- 
calize on the same nucleus. Every 115 as, this ionic state 
is transiently populated by recollision, and in this con- 
figuration the two electrons can easily escape by tunnel 
ionization. We therefore observe a very fast succession 
of Recollision-induced Field-assisted Double Ionization 
(RFDI, see the four double ionization arrows in Fig.[3lc)). 

As mentioned previously, the last efficient pathway to 
double ionization takes place after the end of the pulse, 
and therefore results from an autoionization process. 
This interpretation is confirmed by the snapshot of the 
electronic wave packet shown in Fig.[3](d). For t > 3.0 fs, 
the two electrons escape in opposite directions since the 
ionization is now mainly driven by the electron repulsion. 
This ionization pathway, as well as the RFDI pathway of 
Fig.OJc), are suppressed when we introduce absorbing 
boundaries for the rescattering flux. The origin of both 
ionization channels is therefore to be found in the elec- 
tron recollision. We therefore name the post-pulse ion- 
ization process "Recollision-induced Auto-Double Ioniza- 
tion" , or RADI. Additional calculations performed with 
fixed nuclei show that RADI also takes place when the 
vibrational dynamics is frozen. This process is there- 
fore not induced by a vibrational coupling [17[ , but is 
certainly due to transiently populated dressed molecular 
states lying above the field-free double ionization thresh- 
old. When the field is suddenly switched off, these highly 
excited states decay rapidly by double ionization. This 
interpretation is confirmed by the fact that RADI disap- 
pears for pulse durations longer than 4 optical cycles. 

Finally Fig.[]Jb) reports the contribution of these last 
two ionization channels to the proton kinetic energy spec- 
trum. Since autoionization occurs at larger internuclear 
distances, the corresponding proton spectrum (dashed 
line) appears at lower energies than the RFDI proton 
peak (dotted line). The total spectrum, shown as a solid 
line with circles, includes the four double ionization path- 
ways of Fig.[2ja). However, the direct and sequential 
contributions ® and ® remain very weak in comparison 
with the RFDI and RADI mechanisms. The asymmetry 
of the fragment kinetic energy release spectrum therefore 
reveals the presence of these last two pathways, induced 
by the rescattering process. Additional calculations show 
that the competition between RFDI and RADI is mainly 
controlled by the pulse duration and intensity. RFDI is 
favored with longer pulses and / or higher intensities. 

To conclude, we have presented a theoretical model 
for the study of the electronic and nuclear dynamics of 
diatomic molecules in intense laser fields. This model 
can be adapted to a particular molecule by a simple ad- 



justment of two softening functions a(R) and /3(R). We 
have applied it to the simplest molecular target H2. We 
have shown that at moderate intensities molecular double 
ionization is mainly driven by the rescattering process. 
This process induces a specific transient ionic excitation 
of the molecule which leads to Recollision-induced Field- 
assisted Double Ionization (RFDI) and to Recollision- 
induced Auto-Double Ionization (RADI). These two dou- 
ble ionization mechanisms initiate the Coulomb explosion 
of the molecule at different internuclear distances, an ef- 
fect which can be observed in the fragment kinetic en- 
ergy release spectrum. These mechanisms could also be 
investigated in the time domain using attosecond EUV 
pulses. Indeed, a time-delayed attosecond pulse may in- 
duce an efficient single ionization of H2 in the energy 
range hco ~ 20 eV, thus washing out the electron correla- 
tion in a few tens of attoseconds. Controlling the time- 
delay between the IR femtosecond pump and the EUV 
attosecond probe should therefore make it possible to re- 
veal experimentally these two mechanisms by removing 
them one after the other. 
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